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The effects of the spin-orbit component of the particle-hole interaction on nuclear response func- 
tions and neutrino mean free path are examined. A complete treatment of the full Skyrme interaction 
\^ • in the case of symmetric nuclear matter and pure neutron matter is given. Numerical results for 
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INTRODUCTION 



The mean field theory of nuclear systems is a well developed tool for the microscopic study of finite nuclei as well as 
infinite matter. Many effective two-body interactions have been determined with the goal of achieving a self-consistent 
(^ ■ mean field description of nuclear ground states and excited states through the Hartree-Fock (HF) and random phase 
■^ \ approximation (RPA) approaches. The Skyrme-type interactions are among the most frequently used effective forces, 
^^ . and this is due to the fact that they are relatively simple and yet, they can give a fairly accurate description of finite 
^; nuclei datadlil. 

Homogeneous infinite matter is just an idealized object but it can be a very useful testing ground for various theories 

because practical computations are easier to carry out than in finite nuclear systems. It is also a good representation 

of the internal regions of stellar objects such as neutron stars. The matter equation of state and the response of matter 

' fj ■ to various external probes are important physical properties. For instance, neutrino mean free paths can be deduced 

^ ' from nuclear response functionsfs, 4|. Therefore, HF and RPA studies of infinite matter have always been carried out 

P^ along with finite nuclei studies. The formalism for computing RPA response functions in infinite Fermion systems 

^ ' with Skyrme-type interactions has been shown in Ref.jJ| where applications were made for symmetric nuclear matter 

• ^H , while calculations for neutron matter have been performed in Ref.|a|. The two-body spin-orbit component has been 

^^ ' ignored in the previous studies of the nuclear response function. Whereas it is true that the spin-orbit term does not 

H ' contribute to quantities as the equation of state of saturated spin systems, one should consider also situations where 

■ - - I an external operator can induce spin oscillations which can manifest in the response function. The aim of this work is 

to give a complete treatment of the full Skyrme interaction in the case of symmetric nuclear matter and pure neutron 

matter. The spin-orbit component of the particle-hole (p-h) interaction couples the spin channels in the response 

function. However, it will be shown that the effect remains small, even at momentum transfer larger than the Fermi 

momentum, while the neutrino mean free paths are marginally affected. 

The outline of the paper is as follows. In Sec. II we present the general method for calculating RPA response 
functions with the full Skyrme interaction, generalizing the method of Garcia- Redo et al.(j| to the case of the spin- 
orbit interaction. In Sec. Ill the results obtained for response functions in symmetric nuclear matter and pure neutron 
matter, as well as neutrino mean free paths are discussed. Concluding remarks are given in Sec. IV. 

II. FORMALISM 

A. Definitions 

A general two-body interaction in momentum representation depends at most on 4 momenta. Because of momentum 
conservation there are actually 3 independent momenta. For the p-h case we choose these independent variables to be 



the initial (final) momentum ki (k2) of the hole and the external momentum transfer q. This is illustrated by Fig.Q] 
We will denote by a = {S,M]I,Q) the spin and isospin p-h channels with S=0 (1) for the non spin-flip (spin-flip) 
channel, /=0 (1) the isoscalar (isovector) channel, M and Q being the third components of S and /. 
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\/ \ 
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FIG. 1: Direct and exchange parts of the ph interaction. 

Let us consider an infinite nuclear medium at zero temperature and unpolarized both in spin and isospin spaces. 
At mean field level this system is described as an ensemble of independent nucleons moving in a self-consistent mean 
field generated by the starting effective interaction treated in the Hartree-Fock (HF) approximation. The momentum 
dependent HF mean field, or self-energy determines the single-particle spectrum e(k) and the Fermi level eCkp). 

To calculate the response of the medium to an external field it is convenient to introduce the Green's function, 
or retarded p-h propagator G^"''(q, cjjki). From now on we choose the z axis along the direction of q. In the HF 
approximation, the p-h Green's function is the free retarded p-h propagator Q: 

r rn,,w^ /(ki)-/(|ki + q|) 

where the Fermi-Dirac distribution / is defined for a given temperature T and chemical potential fi as /(k) = 
[1 -I- e'-'^'-'^)"^)' -^]~^. The HF Green's function Ghf does not depend on the spin-isospin channel a. To go beyond the 
HF mean field approximation one takes into account the long-range type of correlations by resumming a class of p-h 
diagrams and one obtains the well-known random phase approximation Q- The interaction appearing in the RPA is 
the p-h residual interaction whose matrix element including exchange can be written as: 

vX^"'\qMM) ^ {ci + -kuK\i^)\V\q + -k2,-k2'i»')) ■ (2) 

The RPA correlated Green's function G]^p^ (q, a;, ki) satisfies the Bethe-Salpeter equation: 

G^^Uq,u;,k,)^GHF{q,ioM) + GHFiq,u;M) Y. J ^V^;:^'''\qMM)G^RpV^,^M) • (3) 

(a'} 

Finally, the response function x''°'Hl,^) in the infinite medium is related to the p-h Green's function by: 



Xrpa {(1,^)^9 J 7CT ^RPA (q, ^, ki ) , (4) 

where the spin-isospin degeneracy factor <? is 4 for symmetric nuclear matter and 2 for pure neutron matter. The 
Lindhard function xhf is obtained when the free p-h propagator Ghf is used in Eq. I^J- 

In the following parts of this work we will often deal with integrals similar to those appearing in Eqs. H3I4|I and we 

will adopt the notation (VC^" CIppAJ, {GpPj^) for such quantities. 

B. The p-h interaction 

The central component of the p-h interaction can be written in the general form: 



'(<Z,ki,k2) = Sia,a')lw^''^+wt\ki+kl]~2Wt~>Y^l^^T.^l^''^)^^^('^^)] ' 



where the W^" are combinations of the Skyrme parameters {ti, Xi) and of the transferred momentum q. Their detailed 
expressions are given in Refs.Q and for the symmetric nuclear matter and pure neutron cases, respectively. One 
can note that there is no coupling between the different spin and isospin channels. The general case of matter with 
an arbitrary neutron-to-proton ratio has been studied in Ref.(3|. 

The Skyrme interactions also contain a zero-range spin-orbit term |j|. It has the form iWso{cri +(12) ■ [k' x (5(ri2)k], 
where Ci is the spin operator of particle i, and k and k' are the relative momentum operator of the particles acting to 
the right and left, respectively. To calculate the contribution of this term to the p-h interaction one has to evaluate 
the matrix element (EJ of this spin-orbit interaction. As this term is density-independent there is no rearrangement 
contribution and the result is just adding the following term 

- S{I, /') w{I) J—qW,o\ SiS, l)SiS', 0)Ms kiYi-Ms (h) - kan-Mslfe) 



+S{S,0)S{S',l)M's 



kiYiM'Jh) 



k2yiM' (fe) 



(6) 



to Eq. 0. We have defined w{I) = 2 + (—1)^ in the case of symmetric nuclear matter, and w(/) = 2 for pure neutron 
matter. The effect of the spin-orbit component is to couple both 5" = and 1 channels. 



Response function 



To obtain the RPA response function of Eq. I^J one has to calculate the correlated Green's function G 
technical details are given in the Appendix. The response function can then be written in the form: 
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(7) 



In this expression kp is the Fermi momentum while m* denotes the nucleon effective mass. The functions xo, X2 and 
X4 are generalized free response functions, defined as ,5] 
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(a) 



with Xo = {Ghf), and xhf = 9Xo- 

The coupling between the two spin channels appears implicitly in the function W^"' which can be expressed in 
terms of the quantities (3i introduced in Ref.Q (their definitions are recalled in the appendix). One obtains: 

132-133 



VFi^") = Ty^(") + G'^''^w^{I)Wl q 



.("') 



(9) 



where a' is defined with respect to a as S" = 1 — 5, /' 
W2 does not depend on them), and 



I (the third components M' and Q' are irrelevant here since 



C-(") = 1 if 

1 9 

= -A/2 
2 



S = 
if S: 



(10) 



If we replace W^""' in Eq.Q by W^ we obtain the results of Ref.Q, as it should be. It is worth noticing that 
the spin-orbit interaction induces a complex coupHng between the 5=0 and 5=1 channels. This coupHng is seen in 
Eqs.(Al-A3) of the appendix. The /3i's are complex and therefore, W^ is a complex function of q and w. 

Finally, the quantity of interest is the dynamical structure function S^"'{q,uj). At zero temperature, it is just 
proportional to the imaginary part of the response function at positive energies. At finite temperature T, one can 
also relate it to the response function by the detailed balance relation and obtainQ: 

lImx("'(q,w,T) 
■~ 1 - e-'^/^ ' 



5(")(q,tu,T) 



(11) 



where the energy u can be either positive or negative. 



III. RESULTS 



We apply the above formaUsm to the case of neutron matter. The calculations are performed with the Skyrme 
interaction SLy230b which is designed for reproducing the neutron matter equation of state[l3|- We can make a 
global assessment of our numerical accuracy by calculating the energy-weighted integrals of our dynamical structure 
functions, and compare them with the energy- weighted sum rule obtained by the double commutator method [ll| . 
Indeed, the sum rule value must not depend on the spin-orbit interaction. 



A. The spin-orbit induced interaction 



r(") 



interaction parameter W^ is real and independent of w and of the temperature, the new function W^' is complex 



q = 2 fm 



1 



The function Wl describes the coupling between spin channels induced by the spin-orbit interaction. While the 

iteraction parameter Wi is real and independent of lo ar 
and it depends on both w and T through the /32, /?3 terms. 
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FIG. 2: The interaction parameter Wj (q,a;) as a function of lo, for two different values of the momentum transfer. Solid and 
dashed lines correspond to S = and S = 1 channels, respectively. The results are for neutron matter at T = and p=po- 

In Fig. 12 are plotted the real and imaginary parts of Wf" (q, w,T), as a function of w for two different values of q. 
The density has been fixed to the value po of saturation density of symmetric nuclear matter, and the temperature is 
T = 0. It is interesting to notice that the w-dependence of the difference (32 — Pz has the same symmetry properties 
as the response function, namely, the real (imaginary) part is symmetric (antisymmetric) with respect to w=0 for a 
fixed value of q. This property is fulfilled in spite of the fact that each of the functions (32 and (3^ does not satisfy it 
separately. Such symmetry properties are seen on Fig. [2 



When uj -^ ±oo, both (3i functions of Eq.JHJ go to zero and W^' tends to Wi . Therefore, the amplitude of the 



M) 



oscillations of the curves in Fig. |2 show the deviations of Wi with respect to Wf" . As for the real parts, these 
deviations are most visible around uj=0. For the imaginary parts, they are always zero at uj=0 and the deviations 
become significant at larger values of \uj\. In any case, one can see that the change of the p-h interaction due to the 
spin-orbit force seems to be relatively small, except at higher values of the transferred momentum q. This increase 
with q refiects the q* power explicitly appearing in Eq.ljHJ. 



B. Dynamical structure functions 

We now turn to the effect of the spin-orbit interaction on the dynamical structure function S"*^"' (q, uo). In Figs. l3l4l are 
plotted the values of S^"'\(\,uj) as functions of w, calculated at zero temperature and at T=20 MeV. The calculations 
are done at the fixed values of momentum transfer q and neutron matter densities indicated in the figures. We show 
separately the S=Q and S=l cases, calculated with and without the spin-orbit force. A general observation is that 
the effect of spin-orbit interaction increases as q increases. This can be easily understood according to the previous 
analysis of the interaction parameter W^' \q,uj). 

We first discuss the T=0 results of Fig. El For the smaller value of q the effect increases slightly with increasing 
density, this effect being more visible in the S=Q channel. For the larger value of q the increasing effect with increasing 
density becomes more dramatic in the S=l channel. The reason for this behaviour is due to the fact that, ai p — po 
and g=2 fm~^ we are approaching the point of instability in the S=l channel. This instability is characteristic of 
Skyrme-type interactions [l3|. Indeed, the response function depends not only on W^ but also on W2 as shown 



in Eq.lO. Now, the interaction SLy230b gives W2^-°=163.55 MeV.fm^, W^ 



=1 



-163.55 MeV.fm^. From Fig.Elone 



also sees that the spin-orbit force contributes some important negative amount in the S=l channel, thus enhancing 
the tendency to instability. One can conclude that the effect of the spin-orbit force is generally small but it can 
become dramatic when one approaches the instability region of the S=l channel. 
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FIG. 3: The dynamical structure function ^^"'(q, tj) at r=0 as a function of w, for two values of the momentum transfer q 
and at densities p = po/2 and po . Solid and dashed lines correspond to S = and 5=1 channels, respectively. The thin lines 
represent the structure function without spin-orbit interaction. 



The above conclusions remain valid at finite temperature but they are further amplified, as one can see from Fig.^ 
in the case of T=20 MeV. As expected, noticeable modifications of the dynamical structure function are found for 
u) around zero and q larger than kj? . One must keep in mind that the effective interaction parameter W^ (q, uj) 
containing the contribution of the spin-orbit force is temperature dependent through the /32, /Sa functions. In the 
T=20 MeV case the vicinity of the instability point is clearly seen for p — po and q=2 fm^^. 
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FIG. 4: Same as FigEl for T =20 MeV. 
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C. The neutrino mean free paths 

We examine now the effect of the spin-orbit interaction on neutrino mean free paths in neutron matter under various 
conditions of density and temperature. The scattering of neutrinos on neutrons is mediated by the neutral current of 
the electro-weak interaction. In the non-relativistic limit and in the case of non-degenerate neutrinos, the mean free 
path A of a neutrino with initial momentum ki is given by 0, |j| 



l/A(ki,T) = 



Crp 

327r3 



dk. 



M^ 



cos 



0)5(°)(g,T) 



4(3 



cos 



0)5«(g,T) 



(12) 



where T is the temperature, Gp is the Fermi constant, cy {ca) the vector (axial) coupling constant. The final neutrino 
momentum is ka, the four- vector q = ki — k^ stands for the transferred energy- momentum, and cos 6 = ki k^. The 
dynamical structure factors S^^'{q,T) describe the response of neutron matter to excitations induced by neutrinos, 
and they contain the relevant information on the medium. The vector (axial) part of the neutral current gives rise to 
density (spin-density) fluctuations, corresponding to the S = {S = 1) spin channel. 

We have calculated the neutrino mean free path at different densities (po and 2po) and temperatures (10, 20 and 
30 MeV). The energy of the incoming neutrino is chosen to be Ei, — 3T. Results are shown in Table |l| The first line 
(Xhf) shows the results for the neutrino mean free path calculated at the self-consistent mean field approximation, i.e., 
without RPA correlation effects. The spin-orbit interaction has no effect in this case, since we are in a homogeneous 
medium and the spin-orbit force does not contribute to the HF properties. Next, we show the results of the complete 
calculation with or without spin-orbit interaction. For T = 10 MeV, the spin-orbit interaction modify the mean free 
path by only 1%, 3-4% for T = 20 MeV and 5-10% for T = 30 MeV. We thus conclude that the effects of the spin-orbit 
interaction on the neutrino mean free paths are at the level of a few percent. 



IV. CONCLUDING REMARKS 



We have investigated the effects of the spin-orbit component of the p-h interaction Vph on the RPA nuclear response 
functions and their possible consequences on the neutrino mean free paths. This study is carried out in the framework 
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TABLE L Neutrino mean free paths (in meters) in neutron matter, calculated in HF and RPA schemes, 
correspond to the different densities and temperatures as indicated. The neutrino energies are Ei, — 3T. 



The calculations 



of a Skyrme-type, zero-range effective interaction. While the central component of Vph keeps the S=0 and S=l spin 
channels separated, the spin-orbit component couples these channels together. However, within the specific form of 
Skyrme-type interactions this coupling appears only implicitly through a modified interaction parameter Wi (q, uj, T), 
and the calculation of the response function is formally identical to the case without spin-orbit interaction. 

The modified interaction parameter W^ is shown to be complex and it depends on the energy-momentum transfer 
(tjjq) and temperature T. Its behaviour at large uj shows that the effect of the spin-orbit force tends to zero for 
increasing uj. The overall effects on the response functions remain small in neutron matter at densities up to po- 
However, in the specific example of the SLy230b force that we have considered, a pole in the response function and 
hence an instability occur in the S=l channel at p ~ po and q ~ 2 fm~^. In this case, even a small modification 
brought about by the spin-orbit force produces a large change of the S=l response function near the pole. 

As for the T-dependence of the spin-orbit effects, all the remarks made above remain true with increasing T, the 
only difference being that the effects are amplified at higher temperature. Finally, the neutrino mean free paths in 
neutron matter are very moderately affected by the spin-orbit component of the p-h interaction. 

The numerical applications have been presented here for the case of neutron matter, but similar results are obtained 
in symmetric and asymmetric nuclear matter. We also note that the zero-range nature of the spin-orbit force studied 
here is reflected in the q'' dependence of the modifled p-h interaction, and therefore the fact that the spin-orbit effects 
increase with increasing q would be altered for a finite range interaction. 

Acknowledgments 

This work is supported in part by the grant FIS2004-0912 (MEC, Spain), the IN2P3(France)-CICYT(Spain) ex- 
change program and the Financement Theorie IN2P3 fund. 

APPENDIX A: THE COUPLED INTEGRAL EQUATIONS FOR THE G*,"^^ GREEN'S FUNCTIONS 



In this Appendix we shall show the way to transform the Bethe-Salpeter equation ^ into a set of coupled algebraic 
equations for three integrated quantities depending on G^pj^. To alleviate the notation, we will omit the (q, tj)- 
dependence, and specify only the spin variables, as there is no isospin coupling. The expressions are formally valid 
for both symmetric nuclear matter and neutron matter. The only differences will be in the value of the spin-isospin 
degeneracy factor g of Eq. Q, and the factor w{I) of Eq. 0. 

One must note that the multipole expansion of Ghf (see Eq.Q only involves terms of the type 1l,o(^i)- Therefore 
integrals of the type {f{k)YL^MGHF) or {f(k)YL,MYL',M'GHF) vanish unless M ^ or M + M' ^ 0, respectively. 
This will simplify the response function equations. Other integrals involving Ghf are also needed and they can be 
expressed in terms of the quantities Pi introduced in J3|: 



{Ghf) — /3o , 
(kyi,oGHF) = q 



tp^ ' 



{\^Ghf) = q'/32 , 
{k^Y,,oGHF)^q'^P4 

2 3 



(k^GffF) = q4/35 



{k^\Yi.o\^GHF) = q^i:P3 , (k2|yi,i|2GHF) = q^^ (/32 - Ps) 



Let first consider the S = channel. With the p-h interaction given by Eqs. H5lfi|l . the Bethe-Salpeter equation is 



written as 
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Integrating over ki we get 



.(0), 



(G'n'pA) = Po + wr(3o{GZA) + wr^'P2{G^n'pA) + wrPo{k'G^n'pA)-^W^' 

+ 2^WsoqPo Y. ^^S i^yiM'sG^RPA^) ■ 



-qpi{kY,oG'^>,^) 



(0) 



(A2) 



One can see that the quantity (Gj^pa) we are interested in is coupled to (k'^G^p^), (kYioGpp^), and {kYiM> G}jp^"^ ). 

Two new equations are obtained multiplying Eq. |IA1|I with k^ and kiYio(fci), and integrating over ki. These factors 

are such that there is no contribution from the term {Gp p/ ) entering Eq. |IA1|I . The coupHng with the S* = 1 
channel is thus contained in the last term entering Eq. IJA2II . From the Bethe-Salpeter equation for the 5 = 1 channel 
the following expression is obtained 



|:(MS«^„,o.-.,..y±..,3^-^^|^^^A_ 



\Grpa) 



(A3) 



This means that the effect of the spin-coupling can be simply absorbed in an effective W[ coefficient. The equations 

for S* = 1 channel are obtained proceeding along similar lines. The expHcit expression of Wi is given in Eq.JH}. 
Finally, the system of algebraic equations can be written in a compact form for both channels as 



r(a)„2, 



,2^(a) 



(i - wi'^^p, - wr^'P2) {G'^u) ~ wrM^'G^pi^^) + 2w^r^q/3ii 
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ql3,{k^Gfp^) + (l + 2M/i"^q2/33) J^{kY^^{k)d^l^) = q/3: 
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